Energy Dependent Tunneling in a Quantum Dot 
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We present measurements of the rates for an electron to tunnel on and off a quantum dot, obtained 
using a quantum point contact charge sensor. The tunnel rates show exponential dependence on 
drain-source bias and plunger gate voltages. The tunneling process is shown to be elastic, and a 
model describing tunneling in terms of the dot energy relative to the height of the tunnel barrier 
quantitatively describes the measurements. 
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Gate controlled quantum dots have been used to study 
a wide variety ofphysical phenomena, from correlated- 
electron physics which becomes important when the 
coupling of the dot to its leads is strong, to the coher- 
ence of electron charge and spin states [2j , which can be 
maintained only when the coupling to the leads is weak. 
Observation of these phenomena is made possible by in- 
situ control of the rate of tunneling T, coupling the dot to 
its leads: T can be adjusted over more than ten orders of 
magnitude by changing the voltages applied to the gates 
that define the tunnel barriers of the quantum dot. Re- 
cently, integrated charge sensors [j| have made possible 
a variety of new investigations of T, probing the effect 
of excited states 0, 0, Ej , electron number Tj barrier 
symmetry 0,0] and Coulomb interactions 

In this Letter, we report exponential sensitivity of the 
rates for electrons tunneling on (r o „) and off (r o //) a lat- 
eral GaAs quantum dot as a function of the drain-source 
bias Vds and plunger gate voltage V g . The tunnel rates 
are obtained by performing time resolved measurements 
of a quantum poin t contact (QPC) charge sensor adja- 
cent to the dot and utilizing pulsed-gate tech- 
niques 0, Q . We show that the tunneling is dominated 
by elastic processes, and that the measured exponential 
dependence of T on Vds and V g is in excellent quantitative 
agreement with a model describing tunneling in terms of 
the dot energy relative to the height of the tunnel barrier. 

The device used in this work is fabricated from a Al- 
GaAs/GaAs heterostructure grown by molecular beam 
epitaxy. The two dimensional electron gas (2DEG) 
formed at the AlGaAs/GaAs interface 110 nm below the 
surface of the heterostructure has a mobility of 6.4 x 10 5 
cm 2 /Vs and density of 2.2 x 10 11 cm -2 [lj|. An electron 
micrograph of the gate geometry is shown in Fig.^a) and 
is similar to that of Ciorga et al. 18]. Applying negative 
voltages to the labeled gates forms a single dot containing 
one electron coupled to the surrounding 2DEG through 
two tunnel barriers: One between SGI and OG (bl, con- 
necting the dot to lead 1) and the other between SG2 
and OG (b2, connecting the dot to lead 2). Lead 2 is 
grounded, and a voltage Vds is applied to lead 1. Ap- 
proximately the same DC voltage V g is applied to the 
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FIG. 1: (a) Electron micrograph of the gate geometry and 
schematic of the measurement circuit. Unlabeled gates are 
grounded. We measure the resistance of the QPC by sourcing 
a current through the QPC and monitoring the voltage V qpc 
across the QPC. We have verified that the finite bandwidth of 
our measurement flEt does not significantly affect the results 
presented in this work using computer simulations, (b) When 
a voltage bias Vds is applied across the quantum dot a small 
current flows and the charge on the dot fluctuates between 
and 1 as shown in (c). We measure the time intervals t on 
(t ff) that the electron is on (off) the dot using an automated 
triggering and acquisition system described in Ref. (lt^. The 
offset in the trace is caused by the AC coupling of the voltage 
preamplifier, (d) Histogram of t on times from data such as 
in (c). Fitting this histogram yields F // as described in the 
text. 

three plunger gates LP1, PL, and LP2. We have verified 
that a single, not double, dot is formed. We measure 
the dot in a dilution refrigerator with a lowest electron 
temperature of 120 mK. 

Applying a negative voltage to the gate QG2, we form 
a QPC between the gates QG2 and SG2. The resistance 
of the QPC is sensitive to the charge on the dot 0, al- 
lowing us to measure the number of electrons on the dot 
in real time [ill H2I Hij ] . If we apply a drain-source bias 
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Vds across the dot (Fig.^b)), we observe the number of 
electrons on the dot fluctuate between and 1: A typ- 
ical trace is shown in Fig. (He). To measure ^ ff, we 
histogram the times t on that the electron spends on the 
dot (Fig. QId)) and fit to an exponential Ae" r °« to " [3. 
We obtain T on from the time intervals t ff in the same 
manner. Using this technique we measure ^off and T on 
as a function of Vds (Fig-Et a )) for the case in which there 
is either or 1 electrons on the dot at the position shown 
in Fig.E^b). As Vds is made more negative, r o // is seen 
to increase exponentially. 

To model the tunneling rates, we linearize the exponen- 
tial in the WKB formula for tunneling through a barrier 
|Tg| for a small perturbation SE to the dot energy E and 
a small deviation SU of the tunnel barrier potential U(x): 
r ~ Tq exp [— (SU — SE) k], where k and Tq depend on 
the details of the barrier potential. In a simple capacitor 
model for the dot we assume a linear dependence of SE 
on small perturbations 8V g and SVds about arbitrary but 
fixed Vds and V g values: SE — —eadsE^Vds — ea g ESV g , 
where ctdsE is the ratio of the drain-source capacitance 
to the total dot capacitance, and a g E is the capacitance 
ratio for the three plunger gates. Similarly, SU varies lin- 
early as SU = -eadsi/SVds - ea gU SV g , where a ds u and 
a g u give the coupling of Vds and V g to the barrier poten- 
tial. There will, of course, be different parameters (XdsUi 
a g u, and k for the two barriers. Note that T depends 
exponentially on SU — SE, and therefore depends expo- 
nentially on Vds and V g : One can show that this holds 
independent of the particular shape U(x) of the barrier 
potential, or the shape of the perturbation to the po- 
tential induced by the change 5V g or SVds- Using this 
linearization, we can write down equations for the Vds 
dependence of r o /^ and F orl , including Fermi statistics 
in the leads, assuming elastic tunneling, and considering 
only a single orbital state of the dot: 

r„ ff = r 2 , e-^(i-f 2 (E)) (i) 

+r 1 ,„e* ffV *(l -h(E)) 
r o „= v^oe- p2SVd \f 2 (E) (2) 
+r ] ri,oe l3l5V *° h(E) 

Here /3i )2 = Ki,2\adsU ia - a ds E\, E = -ea ds EVds - 
ea g E^V g is the energy of the one-electron state rela- 
tive to the Fermi energy Ep, and AV g is the shift in 
V g from the to 1 electron transition. / 2 and /i are 
the Fermi functions of the two leads /2(E) = f(E), 
fi(E) = f(E + eVds), and rj is the ratio between the 
tunnel rate onto and off of the dot for a given lead when 
the one-electron state is aligned with the Fermi level in 
that lead. We expect that rj = 2 because of spin de- 
generacy IHJ > and use this value in the calculations 
below. 

To understand Fig. Eta), we note that whether T in- 
creases or decreases with Vds depends on whether lead 
1 is better coupled to the barrier or the dot, that is, 
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FIG. 2: (a). F on and T ff as a function of Vds for large neg- 
ative Vds- The solid line in the upper panel is based on a 
theoretical fit to the data discussed in the text. (b). Dif- 
ferential conductance vs. Vds and V g , showing the to 1 
electron transition. The tunnel rates for this case are made 
large enough so that the differential conductance can be mea- 
sured using standard transport techniques. The data shown 
in (a) are taken at the position of the dashed line. 

whether etdsU or atd S E is larger. Since bl is closer to lead 
1 than the dot, and b2 is farther from lead 1 than the 
dot, it follows that atdaUi > &dsE > a ds u 2 (see Fig-Hfb)) 
[2l| . Therefore, tunneling through bl (b2) increases (de- 
creases) exponentially with increasing Vds- This is re- 
flected in the signs of the exponentials appearing in Eq.0 
and Eq. 

The solid line through the the r o // data in Fig. Eta) 
is a fit to Eq. which for large negative bias reduces to 
F2,oe~^ 2SVds ■ The rate r o „ is shown as a function of Vd s 
in the lower panel of Fig. Eta). At the two points marked 
"e" in the figure T on increases rapidly as the Fermi energy 
in lead 1 is aligned with an excited orbital state of the dot 
0, giving the 1st and 2nd excited state energies to be 1.9 
and 2.9 meV, respectively, above the ground state. The 
higher-energy states are better coupled to the leads and 
thus r o „ rises rapidly when these states become available. 
These energies can also be measured, with larger tunnel- 
ing rates through bl and b2, using standard transport 
techniques (Fig.Etb)), and the results are consistent. 
We note that r o // does not show any special features at 
these points: This is because the electron decays rapidly 
out of the excited orbital states and subsequently tunnels 
off the dot from the ground state [|| . We therefore con- 
tinue to use Eq. 2] when there are multiple orbital states 
in the transport window. 

In the regions between the points marked "e" , T on is 
seen to decrease exponentially as Vd s is made more neg- 
ative, as expected from Eq.El Note that this decrease in 
r o „, with increasingly negative Vds, occurs even though 
the number of electrons that could tunnel onto the dot 
inelastically from lead 1 is growing, presenting strong ev- 
idence that the tunneling is predominantly elastic, domi- 
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nated by states very close to the dot energy E. There is, 
however, an apparent flattening of T on above the extrapo- 
lated exponential decrease near Vd s = -4 mV, close to the 
second excited orbital state. We find that this line shape 
is consistent with broadening of the second excited state 
by a lorentzian of full- width at half-maximum 7 ~ 10 
/ieV. Calculated line shapes are shown for the broaden- 
ing of the first excited state in Fig. 01 and are discussed 
below. 

If a square barrier is assumed, one can compute 
an effective barrier height U% and width wi for b2 
from the fit in Fig. EJa). For a square barrier T = 
j oe -2w 2V /2 m * (u 2 -E)/h* []]| Linearizing the square root 
in the exponential and estimating otdsE ~ d-dsUz ~ &dsE 
and fo ~ Ei s /h ~ 1 THz, where Ei s is the level spacing 
of the dot, we obtain W2 ~ 130 nm and U2 — Ep rs 5 
meV at Vds = Vq. These values are only logarithmically 
sensitive to fo and thus depend very little on our esti- 
mate of this parameter. The voltages we apply to the 
gates to create the tunnel barriers are the same order of 
magnitude as the voltages at which the 2DEG depletes 
and thus it is reasonable that U2 — Ep is found to be 
comparable to the Fermi energy (Ep w 7.7 meV). The 
value for w 2 is also reasonable given the dimensions of 
our gate pattern and heterostructure. 

We next examine, in Fig. the dependence of T on 
and T ff on Vds for both positive and negative Vds- The 
data are taken with V g adjusted so that E is ~ kT away 
from the to 1 electron transition at Vds = 0. The solid 
lines in Fig. [3] are calculated from Eq. ^ an d Eq. |3 and 
are in good agreement with the data. We note that the 
value of f}\ is very close to the value of (3%. If the height 
and width of bl and b2 are comparable one can show 
that B\ ~ adaCJl — "^^-82, and it is therefore expected that 

At Vds ~ —3 mV in Fig. [21 the first excited orbital state 
is aligned with the Fermi energy in lead 1. For slightly 
less negative values of Vd s , where the excited state is just 
above the Fermi energy, T on appears enhanced above the 
solid line calculated from Eq.|3 The dashed line in Fig. El 
includes broadening of the first excited orbital state by 
a lorentzian of full- width at half- maximum 7 = 13 fJ-eV, 
corresponding to a lifetime broadening of r e i = 50 ± 30 
ps. Emission of acoustic phonons can lead to excited 
state lifetimes ~ 100 ps, but for our device parameters we 
exp ect much slower relaxation from this mechanism [Til 
|23j,|2J|. T on also deviates from the solid curve for Vds It 2 
mV: This deviation may be caused by broadening of the 
first excited state as well. These deviations could also 
be caused by inelastic processes, which might begin to 
contribute significantly to T on when the electron energy 
is sufficiently far below the Fermi energy. 

We now turn to the dependence of the tunnel rates on 
the plunger gate voltage V g (see Fig.0J). For these mea- 
surements, the barriers are tuned so that the tunneling 
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FIG. 3: r on (closed squares) and T a ff (open circles) as a 
function of Vds ■ The solid and dashed curves are calculations 
described in the text. The step features near Vds = result 
from the Fermi distribution. 

through bl is negligible. In the region near AV g = 
(where E ~ kT), the electron hops on and off the dot 
spontaneously because there are both electrons and holes 
in the lead at these energies (Fig. IHb)). In this region 
we measure the rates r on and T ff in the same way as 
described above. In the inset to Fig. Ufa) we compare 
the probabilities that the electron is on and off the dot 
Pan = T on /(T on + T of f) andp ff = T of f / (T on + T of f) to 
f( E , V) = i + leE/kT and 1 ~ f(E, v), respectively, with T 
= 120 mK [ff. 

A pulsed technique (see Fig. ^c)) is used to measure 
T a ff when AV g is made sufficiently negative that ther- 
mally assisted electron tunneling ceases. We begin with 
the electron energy, determined by AVg, well above Ep. 
We then apply a positive voltage pulse to the gate LP2 
to bring the one-electron energy level near Ep, so that 
an electron can hop onto the dot. A short time after the 
pulse the electron will hop off the dot because it is above 
Ep. We record the time at which this occurs, At, mea- 
sured relative to the end of the gate pulse. This process 
is repeated and we make a histogram of the number of 
tunnel-off events vs. At; an exponential fit to this his- 
togram yields T Q ff. An analogous technique is used to 
measure T on when the one-electron state is well below 
Ep. 

To understand the V g dependence of T, we note that 
5U2 = a g u 2 SV g , and a g u 2 < a g E because the three 
plunger gates are closer to the dot than they are to b2. 
The dependence of the tunnel rate through b2 on V g is 
thus very similar to its dependence on Vds (Fig-EJ- Start- 
ing at the far left of Fig.^Ja), we see that T a ff decreases 
exponentially as AV g is made less negative: This is anal- 
ogous to the r o // data shown in the top panel of Fig.E^a) 
and Fig. (for Vds < 0, where the electron tunnels off 
through b2). T Q ff decreases rapidly at the to 1 elec- 
tron transition as the ground state is brought below Ep. 
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FIG. 4: (a) r on and T ff as a function of AV g . Closed 
(open) circles are r a // (r on ) measured by observing spon- 
taneous hopping caused by thermal broadening in the leads 
as depicted in (b). Closed (open) triangles are r o // (r n) 
measured using a pulsed gate technique. Solid lines are cal- 
culations as described in the text, (inset) p on (triangles) and 
Poff (circles) compared to f(E,rj) and 1 — f(E, if) (solid lines) 
respectively as described in the text, (b) Sketch of sponta- 
neous hopping caused by thermal broadening in the leads (c) 
Pulsed technique used to measure T a ff. The top panel shows 
the pulsed modulation of the one-electron state energy E. 
The bottom panel shows a sample time trace. The dashed 
vertical lines indicate when the gate pulse begins and ends: 
The QPC responds to the gate pulse because of direct ca- 
pacitive coupling to LP2. When the electron energy level is 
brought near Ef an electron tunnels on the device (indicated 
by a 1). When the electron level is brought back above Ef 
the electron tunnels off the device (indicated by a 0). 

Continuing into the 1 electron valley we see that r o „ de- 
creases as AV g is made more positive: This is analogous 
to the T on data shown in Fig. (for Vd s > where the 
electron tunnels onto the dot through b2). This decrease 
in T on is further evidence that the tunneling is elastic be- 
cause as the one-electron state is brought farther in en- 
ergy below Ep there are more electrons that could tunnel 
onto the dot inelastically at energies closer to the top of 
the barrier. Y on , however, decreases because elastic tun- 
neling onto the dot happens at a lower energy relative to 
the barrier height. 

We model the data shown in Fig. 0{a) using the same 
equations as for the drain-source dependence (Eq. ^ an d 
Eq. [21 except with T^o set to zero and 5Vd s replaced by 
AVg in the exponents. The solid lines in Fig. 01a) are 



calculated from Eq. Q]and Eq. |2Jwith these adjustments 
and describe the data well. Here again we take r/ = 2, 
although the fit is improved with a smaller value for r\ 
[Io| . We note that the value for /?2 obtained from the AV g 
dependence is smaller than the value obtained from the 
SVds dependence: This is expected because a g E < ctdsE- 

Summarizing, we find that the dependence of T on and 
T a f f on Vds and V g is well described by elastic tunneling 
at a rate set by the difference between the electron energy 
and the barrier height. 
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